Analysis of the dynamical cluster approximation for the Hubbard model 
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We examine a central approximation of the recently intro- 
duced Dynamical Cluster Approximation (DCA) by example 
of the Hubbard model. By both analytical and numerical 
means we study non-compact and compact contributions to 
the thermodynamic potential. We show that approximating 
non-compact diagrams by their cluster analogs results in a 
larger systematic error as compared to the compact diagrams. 
Consequently, only the compact contributions should be taken 
from the cluster, whereas non-compact graphs should be in- 
ferred from the appropriate Dyson equation. The distinction 
between non-compact and compact diagrams persists even in 
the limit of infinite dimensions. Non-local corrections beyond 
the DCA exist for the non-compact diagrams, whereas they 
vanish for compact diagrams. 

Introduction Strongly correlated electron systems are 
often characterized by short range dynamical fluctua- 
tions. Consequently, local approximations like the dy- 
namical mean field approximation (DMFA) success- 
fully describe many of the qualitative properties. How- 
ever, in low-dimensional systems, spatial correlations be- 
come increasingly important and are thought to be re- 
sponsible for e.g. non-Fermi-liquid behavior and d-wave 
pairing in the cuprate superconductors. 
The Dynamical Cluster Approximation (DCA) was intro- 
duced as a technique to include such non-local corrections 
to DMFA This is accomplished by mapping the 

lattice problem onto that of a self-consistently embedded 
cluster, with periodic boundary conditions. The DCA 
may also be viewed diagrammatically as an approxima- 
tion which systematically restores momentum conserva- 
tion at the internal vertices of many body Feynman dia- 
grams, which is relinquished in the local DMFA. Here, we 
investigate one of the central approximations of the DCA: 
that compact (skeletal) contributions to the thermody- 
namic potential are well approximated by their cluster 
counterparts; whereas non-compact (non-skeletal) con- 
tributions should not be approximated by their cluster 
counterparts. Rather, they should be constructed using 
the appropriate Dyson equation. 

General considerations Following Baym a micro- 
scopic theory may be defined by its approximation to the 
generating functional $[G'] defining the thermodynamic 
potential (difference) of the system via 

An^n-no^ -2Ti-[J:G - ln(G/Go)] + $[G] , (1) 



non-compact diagrams are illustrated in Fig. |^. In the 
non-compact diagram, two self energy pieces cr and a' 
are connected with two Green functions. In the compact 
diagram, two vertex parts, F and F', are connected with 
four Green functions. 
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FIG. 1. (a) typical non-compact (non skeletal) and (b) typ- 
ical compact (skeletal) diagrams. 

As shown by Miiller-Hartmann 1^, the DMFA may be 
defined by relinquishing the momentum conservation at 
each internal vertex in <&[G] This conservation is 

described by the Laue function 



A — \ ^ „ix-(ki+k2-k3-k4) _ ATr 



(2) 



In the Dynamical Mean Field Approximation (DMFA), 
momentum conservation is completely abandoned and 
the Laue function Ajjj^ipa = 1- We may then sum freely 
over all the internal momenta labels. 
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FIG. 2. The real lattice cluster (left) and the first Brillouin 
zone (left) divided into subcells. 



where G is the full and Go the bare one particle Green 
function, and S the self energy. $[G] is a sum of all com- 
pact (skeletal) closed connected Feynman diagrams. The 
other contributions to the thermodynamic potential in- 
corporate non-compact diagrams. Typical compact and 



The DCA is constructed to systematically restore the 
momentum conservation at each internal vertex, by map- 
ping the lattice onto a self-consistently embedded cluster 
problem. We have provided a microscopic definition of 
the DCA through its Laue function. However, to clarify 
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the relation between this microscopic definition, and the 
cluster problem, we must first decompose the lattice into 
clusters, and define the corresponding problem in recip- 
rocal space. Here, the real lattice of N sites is tiled by 
N/Nc clusters each composed of Nc = sites where D 
is dimensionality and L the size of clusters (c.f. Fig. || for 
L — D = 2). We label the origin of the clusters by x 
and the Nc intra-cluster sites by X. So for each site in 
the original lattice x = X + x. In the reciprocal space, 
the points x and X form lattices labeled by k and K 
respectively with — • 27r/L and integer . Then 
k = K + k (see Fig. |). 
In the DCA, we first make the separation 
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^ - — '^ki+k2.k3+k4 ^C<5ki+K2,K3+K4 (3) 



and then set ^/A^c'5ki+k2,k3+k4 - ^' ttiat 



^DCA 



c"Ki+K2,K3+K4 1 



(4) 



which indicates that the momentum is partially con- 
served for transfers between the cells. 

In this paper, we consider the approximation made 
through the substitution A AocA in the compact and 
non-compact contributions to the thermodynamic poten- 
tial. Whenever the substitution is made, all internal legs 
are replaced by the coarse grained Green function defined 
by 



N, 



(5) 



The corresponding estimate of the self energy will then 
necessarily be a function of K and (complex) frequency 
z. The DCA estimate of the lattice Green function is 
then given by 



G'(k, z) = 



1 



Ck + M 



^DCA 



(K,z) 



(6) 



It is of importance to note that by using Eq. ^ we have 
already made the approximation E(/c, z) = S£)Cyi(K, z) + 
0[{Aka)^) where Afc is the size of the coarse graining 
cell shown in Fig. || and a is the lattice constant (chosen 
as unity). We also drop the frequency label from this 
point on for simplicity. 

Coarse graining compact vs. non-compact diagrams 
We investigate the additional approximations associated 
with coarse graining in compact and non-compact dia- 
grams. To do this we will consider diagrams with all 
legs coarse-grained, except those explicitly displayed in 
Figs. ^ and ^. Consider the first non trivial correction to 
the coarse grained non compact diagrams as illustrated 
in Fig.| 
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i)[Ar!„,p]^— J2 fT(Ki,K2)a'(Ki,K2) / 



/oo 
deide2G(Ki, ei)G(K2, 62) [/9„cp(ei, Ki; £2, K2) - 
-00 

p„cp(ei,Ki;e2,K2)] 5k2,Ki , (7) 
where 



(ei,Ki;e2,K2) ^ ^ J2 ^^^^ "^Ki+ki) 



N 



ki,k2 



^i^2 -eK2+k2) ^k2,kl 



(8) 



and 



(ei,Ki;e2,K2) = 2^ '^(^i - <^K,+k^) 



Pncp \ 



ki,k2 



'^(£2 -eK2+k2)- 



(9) 
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FIG. 3. First correction by non-compact diagrams. 



In the above derivations we assumed that the self energy 
is k independent and therefore the entire k dependence 
of the Green function is only through the dispersion e 
given by 



D 



=K-|-k 



-2t/(2i^)i/2.^cos(K„ + k„), (10) 



which is just the non interacting dispersion of the Hub- 
bard model Hamiltonian with nearest neighbor hoppings 
and dimensionality D. 
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FIG. 4. First correction by compact diagrams, 5'^'[Af2cp]. 

By the same token, in Fig. ^, for the compact part with 
coarse grained F and F' we have 



Ki,K 
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S^^'^lAflc ] ^ — (5k4 Ki+K2-K3r(Ki, K2, K3, K4) X ^"^^ * finite dimensions, we expand the dispersions 

K i< 1 2 3 1 ' 1 exponentials as functions of k and keep the terms 

K3.K4 up to Afc^. The difference (5^'„cp for the non-compact 

,,,, ^ /""^ /"^ /""^ Z""" , .X diagrams becomes 

r'(Ki,K2,K3,K4) / / / / dcide^de^^de^ (11) ^ 

J— 00 J — 00 J — OQ J — 00 ^2 

G(Ki, ei)G(K2, e2)G(K3, e3)G(K4, £4) [pcp(ei, Ki; £2 <5vl/„,p(si, Ki; S2, K2) « -A/c^ x 

,K2; e3,K3;e4,K4) - pcp(ei, Ki; e2, K2; £3, K3; £4, K4) ] , [?7(Ki + K2) - t](Ki - K2)] S1S2 x 

^here exp [i(si£Ki + S2£k2)] , (15) 

^3 where ?7(K) = 1/D cos(K„) . Reversing the Fourier 

Pcpi^i, Ki; £2, K2; £3, K3; £4, K4) = transform then yields for the difference of pncp and pncp 

^1-^2 2 

X (5p„cp(£i,Ki;£2,K2) « -— Afc X 

Sie. - ^K.+kJ ^k.,k.+k.-k3 - (12) [^(Ki + K2) - V{K, K2)] X 

and ^<5(£i -£Ki)^<5(£2 -£K2)- (16) 

Pcp(£i, Ki; £2, K2; £3, K3; £4, K4) = Eq. |l6| demonstrates that the difference between coarse 

graining and not coarse graining in non-compact dia- 

j;^ ^ H^i ~ ^Ki-i-k;i)'^(^2 ^ ^K2+k2) ^ grams yields an error of order Afc^. After tedious but 

ki,k2 straightforward calculation for the corresponding differ- 

^' " ence between coarse graining and not coarse graining 

'^(^3 ~ ^K3-i-k3)'^(^4 ^ ^K^+kJ ■ (1"^) in compact diagrams we obtain Spcp 0{Ak^). Since 

A/c = 2tt/L, for large clusters this error becomes much 
smaller than the error in Eq. |l|. 
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FIG. 5. The compact contribution of the thermodynamic 
potential versus Ak^ using FLEX method for U = 1.57. 

We now define the Fourier transforms of p and p with 
respect to the energy arguments, e.g. 



p{e,K)e""'d€. 

-00 



(14) 



For the non-compact corrections, Eq. O defines 
*„cp(si, Ki; S2, K2) and 'l'„cp(si, Ki; S2, K2)7 respec- 
tively. In order to calculate the difference between "if 
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FIG. 6. non-compact contributions of the thermodynamic 
potential with and without coarse graining for U = 1.57. 

To illustrate the above point we simulate a two- 
dimensional Hubbard model (with local interaction JJ 
and a nearest-neighbor hopping i = 1) using the Fluc- 
tuation Exchange Approximation (FLEX) [Q. We em- 
ploy an elaborate subtraction scheme to correctly deal 
with the high frequency behavior of the Green functions 
and FLEX potentials Q. In Fig. || and Fig. ^ the com- 
pact and non-compact parts of Af2 have been plotted 
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for the interaction U = 1.57 for various cluster sizes L 
and temperatures T. In Fig. it is readily seen that 
for the compact contribution with coarse grained Green 
functions the variation of Slcp over the entire Afc range 
is roughly 10%. In contrast, as shown in Fig. ^, the dif- 
ference between the non-compact contributions with and 
without coarse graining generates a deviations of over 
100%. Note that at very low temperatures (inset in Fig. ^ 
) the deviation from linearity for T < 0.066 is due to the 
correlation length exceeding the size of cluster L. Then, 
the DC A assumption of replacing S(k) by SDC'yi(K) is 
no longer valid. 

Limit of Infinite Dimensions To make contact with the 
original derivations of the DMFA [Qj^], we explore the 
differences between the compact and non-compact graphs 
in the limit of many spatial dimensions, D. For D —> 
oo, we calculate Ap for both compact and non-compact 
diagrams by simply expanding ^' and ^ but this time 
with respect to 1/D instead of Ak^. Letting D —^ oo we 
can show that 

lim *„cp(si,Ki;s2,K2) « *d^oo x 

D—^oc 

exp(^Ak2[,;(K2 - Ki) - 7/(Ki + K2)]siS2) , (17) 

Expanding the exponential in Eq. |l^ 

lim *„cp(si, Ki; S2, K2) « ^d^oo x 



D^oo 



(1 - |^AA:2[ry(K2 - Ki) - 7y(Ki + K2) 



that the Green function falls off exponentially as G{r) ~ 
e~'^/^= where is a characteristic length depending on di- 
mensionality. It turns out that the compact diagrams fall 
off much faster due to having a larger number of Green 
functions relative to the non-compact ones which again 
explains why we coarse grain the Green functions in only 
the compact part of the free energy difference. The ar- 
guments made here may also be applied to the DMFA 
simply by taking Nc — 1. They also may be extended to 
other cluster approaches, such as the Molecular Coherent 
Potential Approximation |12| and its formal equivalent 
for dynamical systems, the Molecular Cluster Dynamical 
Mean Field |l3|. 

Conclusions We justify one of the central underlying 
approximations of the dynamical cluster approximation. 
We both analytically and numerically demonstrate that 
coarse graining the Green functions in non-compact di- 
agrams results in a larger amount of error compared to 
that incurred in compact diagrams. Consequently, non- 
compact diagrams and their contribution to the thermo- 
dynamic potential are not coarse-grained. The distinc- 
tion between non-compact and compact diagrams per- 
sists even in the limit of infinite dimensions. This con- 
curs with previous work on dynamical mean field theory 
and has implications for other cluster approaches. 

We would like to acknowledge Th. Maier, R. Jamei and 
A. Voigt for very useful discussions. This work was sup- 
ported by NSF grants DMR-0073308. 
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and consequently 



(18) 



lim (5p„cp(ei,Ki;e2,K2) 
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-Afc2 X 



[7;(K2 - Ki) - 77(Ki + K2)] — — 
pz3^oo(ei,Ki;e2,K2) + 0(Afc4). 



(19) 



However, for the difference in the compact diagrams as 
— > 00 we find lim£)_,oo 6^ cp — 0. Thus, for all 
Afc, there are non-local corrections to non-compact di- 
agrams while there are none to compact diagrams as 
D 00. This result is consistent with what Miiller- 
Hartmann has shown for the DMFA 1^. Regardless of 
whether the expansion parameter is Afc or 1/D, there 
is a fundamental difference between compact and non- 
compact diagrams which requires different treatments 
of each. Consequently, only the compact contributions 
should be formed from coarse-grained Green functions, 
whereas non-compact graphs should be inferred from the 
appropriate Dyson equation. 

Finally, even though we directly invoked the nearest 
neighbor hopping dispersion in our algorithm, our ar- 
guments can be generalized to other models assuming 
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